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Some examples in support of the conjecture that the horizon area of a near equilibrium black 
hole is an adiabatic invariant are described. These clarify somewhat the conditions under which the 
conjecture would be true. 



, I. INTRODUCTION AND SUMMARY 

o^ : 

On ' Consider the event horizon area of a black hole. Hawking's area theorem |1| indicates that whenever classical fields 
' obeying the weak energy condition are involved, this physical entity must always grow under external disturbances. Be 
O ^. that as it may, it is possible to exhibit an assortment of cases in which the slow application and relaxation of an external 
D ' disturbance does not necessarily mandate an area increase. This is of course in agreement with the identification of 
C/^ the horizon area and the entropy associated with the black hole [Q. After all, in classical thermodynamics entropy 
' is invariant under slow, reversible changes of a system in thermodynamic equilibrium. This implies a similarity of 
horizon area to an adiabatic invariant in mechanics. What is an adiabatic invariant? 
CO • Consider a Hamiltonian system characterized by a Hamiltonian function H{p, q; A) depending on a parameter 
^ ^ A. Suppose that the parameter changes slowly with time. It turns out that in the limit as the rate of change of 
the parameter approaches 0, there is a remarkable asymptotic phenomenon: two quantities, generally independent, 
become functions of each other. For example, consider the motion of a perfectly elastic rigid ball between perfectly 
elastic walls whose separation / slowly varies. In this case the product vl of the velocity of the ball v and the distance 
between the walls turns out to be an adiabatic invariant, e.g. if we increase very slowly the separation and then very 
00 . slowly decrease it to the original value (slow in comparison with the characteristic oscillation), then at the end of this 
On ■ process the velocity of the ball will be the same as it was at the start. Additionally, it turns out that the ratio of the 
' energy H of the ball to the frequency of the oscillation between the walls uj changes very little under a slow change 
of the parameter, although the energy and frequency themselves may change a lot. Quantities such as this ratio are 
^' ^ ' called adiabatic invariants. A more precise definition of adiabatic invariant is |^ 
{3J[), Definition. The quantity I{p, q; A) is an adiabatic invariant of the system 
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; ^=^9^' '^^a^' H^H{p,q-et). (1.1) 

u : 

5^ if for every k > there is an eo > such that if e < eq and < t < 1/e, then 

\I{p{t), q{t)-et) - /(p(0), g(0); 0)1 < n. (1.2) 

Note that this definition allows for the preservation of the adiabatic invariant with a power-law accuracy (with respect 
to the small parameter describing the speed of frequency change), in comparison with the exponential accuracy 
considered to be generic 

In quantum terms, one may comprehend the adiabatic invariance of H /w by the following arguments. Since it 
is assumed that during the adiabatic change of the system the perturbations imposed have frequencies very small 
compared to w, transitions between states of successive quantum numbers are strongly suppressed. Now, for example, 
for a harmonic oscillator in a stationary state labeled by quantum number n, H /lj ~ [n + l/2)h. Therefore, the ratio 
H/lo is preserved. 

A general proof that the horizon area is an "adiabatic invariant" will enable the utilization of Ehrenfest's principle [|] : 
any classical adiabatic invariant corresponds to a quantum entity with discrete spectrum. Ehrenfest showed that for 
a quasiperiodic system, all Jacobi action integrals of the form / = §pdq are adiabatic invariants. This principle is 
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the backbone of the Bohr-Sommerfeld theory- "old quantum mechanics" , where all Jacobi actions, invariant under 
adiabatic changes, are quantized in integers. 

Even so, application of this principle to the problem at hand or to any other problem does not depend on whether the 
adiabatic invariant can be expressed as a Jacobi action integral or not. In fact, since the Hamilton- Jacobi formalism 
for black holes is not very well developed, it is not even clear that the horizon area of a black hole can be so expressed. 
First attempts in this direction were made in the framework of the membrane paradigm In this framework the 
problem of slow time dependent perturbations of a black hole is the basis for the membrane description of black hole 
evolution. This describes the horizon as a viscous conductive two dimensional surface. Recently Parikh and Wilczek 
have calculated an action for black hole membranes using the Euclidean formulation of Teitelboim ||]. Hence, as 
matters stand we find it useful to seek insights into the problem using a different approach which is based on the 
Newman-Penrose formalism as presented in the next section. 

Accordingly, application of Ehrenfest's principle will then lead to the spectrum for the horizon area operator. On 
that account, the object of this work is to present a collection of examples in which the horizon area of the black 
hole shows signs of being the analogue of a mechanical adiabatic invariant (not necessarily in the form of a Jacobi 
action integral). Each example is decomposed into a static and dynamical problem. As can be easily understood, in 
the static case there cannot be any energy flux through the horizon, and thus no area increase. However, the static 
case serves as a limiting case for the dynamical problem as we take the adiabatic parameter to zero (it is not clear 
beforehand that this limit actually exists). These examples suggest the existence of a theorem which would state 
that, classically, under suitably adiabatic changes of a black hole in equilibrium, the area of its event horizon does 
not change. This would provide a formal motivation for quantizing the black hole in the spirit of the "old quantum 
theory" . This approach has lately received increased attention ^ . 

The paper is organized as follows. In Sec. Q we present the basic theory required to answer the question regarding 
the circumstances under which the black hole horizon area remains unchanged under the influence of an exterior 



perturbation. In Sec. [II we discuss the generalization of the result obtained by Bekenstein pX| ] regarding the influence 
of scalar charges in the Schwarzschild black hole exterior on the horizon's area to the case in which the black hole 
is electrically ch arg ed. The extreme case of the previous example is discussed in Sec. As a further example we 
consider in Sec. M the possibility that the charges are distributed inside a de Sitter universe with a Schwarzschild 
black hole embedded inside it. The influence of the perturbation on the cosmological horizon will be also discussed. 
In Sec. VI we study in brief the possibility that the hole itself is rotating. Finally, in Sec. VII we examine the influence 
of an electromagnetic radiation on the horizon area of a Schwarzschild black hole. 



II. THE BASIC IDEA 



Here we shall lay out the basic theory needed to answer the question regarding the circumstances under which the 
black hole horizon area remains unchanged under the influence of an exterior perturbation. 

Consider a small patch of event horizon area SA; it is formed by null geodesic generators whose tangents vectors 
are P = dx"/dX, where A is an affine parameter along the generators pl| . The complex conjugate null vectors m" 
and to" together with l" and a fourth real vector n" orthogonal to them which satisfies I'^Ua — 1, constitute the 
Newman-Penrose tetrad which lies at the horizon, to" and to" are orthogonal to 1°^ and satisfy mf^rfia = — 1. 

By definition of the convergence of the generators p, 5 A changes at a rate 

^5A = ~2p5A. (2.1) 
uA 

Now the convergence, defined hy p = la]p m^m^, changes at a rate given by the optical analogue of the Raychaudhuri 
equation 12 1^ 



p^ + aa + inTa^n^, (2.2) 

where a = la;pm"m^ is the shear of the generators and Ta/s is the energy momentum tensor of the matter in the 
exterior of the black hole. The shear, which measures deformation of the horizon, evolves according to 

^^2pa + Ca^p^s rmH^m\ (2.3) 

where Cafs-ys is the Weyl conformal tensor . Utilizing the orthogonality relations of the tetrad base and Einstien's 
equations one may calculate the Weyl conformal tensor on the horizon 
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CaMs rmf^Pm^ = Ra0js l°'m>^Prn^ + AttT^^p l^f. (2.4) 



It will be noticed that to keep the horizon area constant requires p — Q which by Eqs. ( |2.2| ) (2.3) implies that both 
Cap-yS l^m^Pm^ and Tap vanish at the horizon. Thus, if initially the convergence and shear were equal to zero, 
then they would remain so. Vanishing of Cap-yS l°'m^Pm^ requires that the geometry be quasistationary to prevent 
gravitational waves, which are quantified by C^^^^, from impinging on the horizon. Thus with a quasistationary 
geometry, preservation of the horizon area requires 

Tap ri'^ = on the horizon. (2.5) 

No te th a t in the dynamical case one should expect a total area increase which can be calculated by first integrating 



Eqs. (2.2) (2.3) and then substituting the result in Eq. ( p.l[ ). However, since we are mainly interested in the adiabatic 
character of the horizon area it is sufficient to prove that the rate of change of the horizon area, given by {d 6A/dX) /6A 
is smaller than the rate of deformation, g iven by cr, (note that the comparison is between two quantities with the 



same dimensions). Now, according to Eq. (2.1) {ddA/dX)/SA is proportional to p. Hence, it is sufficient to prove that 
the differential change in p is smaller than the differential change in a. Evidently the entire evolution of the various 
quantities in the problem is governed by the magnitude of the Riemann tensor and the energy-momentum tensor on 



the horizon. Therefore, according to Eqs. (2.2) (2.3) all that is left to verify in order to prove our claim is that the 



term T^p I'^l^ is of higher order in the adiabatic parameter than the term Rap^s l^m^Pm^ . 

A quick way to see that it is indeed a true claim goes as follows. The distribution of static perturbators around a 
static black hole should not cause an area increase (this is but a simple consequence of staticity). This can be realized 
by noting that the zeroth order term in the adiabatic parameter of Rap-yS l^mf^Pfn^ must vanish. However, in general 
the first order term would not vanish. For example, such terms appear when a contraction is made between the first 
order correction to the Riemann tensor and the zeroth order corrections to the tetrad base or vice versa. On the other 
hand, the contraction of the energy-momentum tensor with the generators of the horizon must be an even function of 
the adiabatic parameter. This claim is based on the demand that this term should be a semi-positive definite quantity 
as indicated by the weak energy condition, which is assumed to be obeyed by most of the classical forms of matter 
found in nature. The lowest order available is of course the second one. 

Contrary to the common belief which necessitates an increase of horizon area when changes in the black hole take 
effect, we shall here examine a variety of processes for which the conditions that keep the horizon area unchanged 
occur naturally. The rule that seems to emerge is that quasistationary changes of the black hole induced by an external 
influence will leave the horizon area preserved. This means an "adiabatic theorem" for black holes must exist. 



III. REISSNER-NORDSTROM BLACK HOLE DISTURBED BY SCALAR CHARGES 

In this section we extend the result obtained by Bekenstein regarding the influence of scalar charges in the 
Schwarzschild black hole exterior on the horizon area to the Reissner-Nordstrom case. As mentioned already in the 
introduction, this example, as for the following examples, is decomposed into a static and dynamical problem. As can 
be easily understood, in the static case there cannot be any energy flux through the horizon, and thus no area increase. 
However the static case serves as a limiting case for the dynamical problem as we take the adiabatic parameter to 
zero. 



A. The Static Problem 



Consider a Reissner-Nordstrom black hole with the metric 

ds^ = -(1 - 2M/r + Q^/r^)dt^ + (1 - 2M/r + Q^jr^Y^dr^ + r^{dQ^ + ^viQ^d^'^). (3.1) 

The horizon is located at = M + \J — Q"^ . If the scalar's sources are weak, one may regard the scalar field, $, 
as a quantity of first order, and conduct a perturbation analysis. The scalar's energy-momentum tensor, 

= Va$V^$ - i(5f V-y^V^^ (3.2) 

will be of second order of smallness. Thus to first order the metric ( |3.lD is preserved. 

To begin with we shall consider first the case in which the scalar charges are static. Then, the field equation outside 
the scalar's sources can be written in the form 
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§-({r'-2Mr + Q^)f)-L^^ = 0, (3.3) 



where is the usual squared angular momentum operator. Following an almost identical derivation as in for 
the Schwarzschild black hole case, one finds the following physical solution 



(- 



e=o m=- 



where the are the familiar spherical harmonic and Pe{x) are the well-known Legendre polynomials |T5[| . The 
(complex) coefficients Cgm permit us to match the solution to every distribution of sources by the usual methods. 

An explanation of the term "physical solution" is in order. For the solution to be physically acceptable, 4> must 
not induce any divergences in the invariants of the curvature via Einstein's equations. This means that the invariant 
Vq,$V"$ must be bounded everywhere (any other invariant is proportional to some power of this one). Hence, if the 



other independent solutions of Eq. (3.3), furnished by the Legendre associated functions Qe{x) (which are singular 
for X = 0) were kept, this condition would be violated. 

Keeping in mind that we are interested in the behavior of the expression Tapl'^l^ near the black hole's horizon, we 
turn now to determine the components of the vector P in the new geometry. Recall that P is the tangent to the 
null generator of the horizon. Furthermore, P is future directed as well as outgoing. Actually, there is no need for 
determination of the components of P. Any 3D"hypersurface of the form {Vt, r — const} has a tangent r" = with 
norm —(1 — 2M/r + Q^/r^) as well as the normal 770, — da{r — const) = Sa"^ with norm (1 — 2M/r + Q^/r^). The 
linear combination of these two vectors, N" = r" + (1 — 2M/r + jr^^rf is obviously null, and as r ^ r>^, both 
its covariant and contravariant forms remain well defined. Hence, it must be proportional to Note that just like 
P, N'^ is future directed as well as outgoing. Applying the new found N'^ and substituting the energy- momentum 



tensor for the scalar field, Eq. (3.2) one is able to write 



lim T^^jA^^iV'' = lim (A^"V„$)2 = lim {{\ - 2M jr ^ lr^)d^ j drf = (3.5) 



Equation (3.5) implies that the horizon area is left unaffected by the presence of the static perturbations in the outer 
region of a Reissner-Nordstrom black hole. This is essential if the area is an adiabatic invariant under dynamical 
perturbations of the black hole as is shown in the next subsection. 



B. The Time Dependent Problem 



The time dependent scalar equation in the background Reissner-Nordstrom spacetime is 



In analogy with Eq. (3^) we now look for a solution of the form 



<J D—r\ ™— D 



)e-"^*. (3.7) 

f=0 m=-l 



In terms of Wheeler's generalized "tortoise" coordinate [Q, defined by dr* jdr = [1 — 2M / r + Q'^ / P)^'^ for which the 
horizon resides at r* — —00, the equation satisfied by the new radial function Hii^{r*) = rfi{uj,r) is a Schrodingcr 
type equation 



dr 

V.,,0^(l-^ + ^)(^-^ + ^<^)^ (3^8, 

Here Vi{r*), with r considered a function of r*, is a curvature induced potential barrier p^ . The exact position of 
the peak depends on M, Q and £ in general. When Q = 0, the peak occurs between r = 3M /8 (^ = 0), and r — 3M 
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{£ oo); when \Q\ = M the maximum is at r = 2M for all i. The height of the potential is inversely proportional 
to M^: for Q = 0, the maximum is Vo^max = 27/(1024M2), Ve^max = ii^ + l)/[27M^) as £ ^ oo; for \Q\ = M, 
Ve^max = (1 + 2£(^ + 1))/(32A/^). Note that as we approach the horizon the potential decays exponentially with a 
vanishing value on the horizon. 



The analogy between Eq. (3.8) and the Schrodinger eigenvalue equation calls for the following analysis of the effects 
of distant scalar sources on the black hole horizon . Waves with "energy" on their way in from a distant source 
run into a positive potential barrier. Therefore, waves with any £ and smaller then the maximal hight of the 
barrier, coming from sources at r 3> 3M, have to tunnel through the potential barrier to get near the horizon. In 
due course, the wave amplitudes that penetrate to the horizon are small fractions of the initial amplitudes, most of 
the waves being reflected back. This means that adiabatic perturbations by distant sources perturb the horizon very 
weakly. Thus one would not expects significant growth of the horizon area from scalar perturbations originating in 
distant sources. 

What if the scalar's sources are moved into the region r-u < r < 3M inside the barrier? They will now be able to 



perturb the horizon; do they change its area? To answer these questions, we look for the solutions of Eq. (3.8) in the 
region near the horizon where the potential is small compared to w^. 

But before undertaking this task one should note that as uj tends to zero the sources must also be moved closer 
and closer to the horizon since the inside of the potential barrier moves in towards the horizon. This suggests that 
there are in fact two time scales in the problem; one due to the time scale of the sources and a second one due to the 
rate of change of their position. However, since the potential decays exponentially as we come closer and closer to the 
horizon this means that we can lower lo without cha ngin g the location of the sources very much. 



Consequently we turn to find the solutions to Eq. (3.8). According to the theory of linear second order differential 
equations they are of the form 

He^{r*) = exp(±icjr*) x [1 + 0(1 - 2M/r + Q^/r^)]. (3.9) 

The Matzner boundary condition that the physical solution be an ingoing wave as appropriate to the absorbing 
character of the horizon selects the sign in the exponent as negative. Hence the typical term in $ is 

1 + 0(1 ~ 2M/r + QV^') p^(^^g ^pg^. ^p = uj(r* +t)-m^. (3.10) 

r 

As can be easily verified, invariants of the curvature associated with this solu tion f all proportional to (VQ<i>V"<i>)'^] 
are bounded. For example, consider a $ composed of a single mode as in Eq. ( 3.10 ). An explicit calculation on the 
Reissner-Nordstrom background using dr* /dr = (1 — 2M/r + Q^/r^)"^ gives 

„ m^Pp^ siv? %h ( dPf \ ^ cos^ tb iv sin(27/') , 

V^ci>V"$(x /. + (-1) + 3.11 

sm 6 \ dO J r-^ 

where " • . . " here and henceforth denote terms that vanish as r goes to ry^. This expression is obviously bounded 
at the horizon. By induction it is possible to show that a $ which is the sum of several modes is bounded on the 
horizon. 

Finally we turn our attention to the question whether the horizon area is left unchanged by the perturbation. As 
with the static case, the expression which governs the extent by which the horizon area increases is proportional to 

TapN"'N'^ = {d<^/dt + (1 - 2M/r + /r^)d<^ / drf (3.12) 
If one now substitutes a made up of a single mode like in Eq. ( 3.10| ), one finds out that 

T.,n^oc "'^;f'^ +.... (3.13) 

In conclusion, on the one hand the amount by which the shape of the horizon is perturbed when the scalar field 
sources are moved inside the potential barrier does not vanish as the perturbation ch anges slo wly. This quantity is 
measured by the complex shear as defined in the previous section. According to Eqs. (|2.2| )(pT^), the evolution of the 
shear is governed by the value of the Weyl tensor on t he h orizon. This can be expressed by the Riemann tensor and 
the energy momentum tensor on the horizon [see Eq. ( [2.4| )]. Now, the metric perturbation, via Einstein's equations, 
in general must be linear in t he m agnitude of the invariant Vq$V"$ (after all Einstein's equations have T^f^ as 



source, not To^T^^). By Eq. ( 3.11 ) this perturbation is of order 0(a;°) generically, and of 0{lju) in the monopole 
case. Clearly, the zeroth order of perturbation 0{uP), corresponds to the static case which was already dealt with 
previously. Therefore, the dynamic deformation (to be compared with static deformation) is of order O(ti^). On the 
other hand we have shown that the rate of change of the horizon area is of order 0{uP'), so that as w ^ the change 
in area becomes entirely negligible on the scale of the perturbation. 
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IV. EXTREME REISSNER-NORDSTROM BLACK HOLE DISTURBED BY SCALAR CHARGES 



It is a well-known fact that the area of an extreme Reissner-Nordstrom black hole is not invariant under the 
absorption of a point particle . But as we shall prove now the horizon area of an extreme Reissner-Nordstrom 
black hole is indeed an adiabatic invariant under scalar field-black hole interaction. 

A. The Static Problem 

The scalar field equation in the extreme Reissner-Nordstrom background for a time independent situation is 
with the general solution 



«=0 m=-i ' / 



(4.2) 



The complex coefhcients Aim and Bi„i enable us to match the solution to any distributions of sources by the usual 
methods. Now, since we are in search of a solution which is regular on the horizon [recall that invariants such as 
(Vq,$V"$)'^' must be well behaved there], we must make all the coefficients Aim vanish. Thus in the inner region of 
the black hole exterior (inwards from the sources) we have, to first order in perturbations theory, the exact solution 

oo e 

'^ = ^E E Bem{r-MYYem{e,'f). (4.3) 

t=0 m=-i 

As with the non-extreme case we shall employ the vector A^". In the new geometry, = t" + (1 — M/rYrj". The 
change in the horizon area as measured by 

lim Taf3N°'N^ = lim (Af"V„$)2 = lim ((1 - M/rfd'^/drf (4.4) 

r^M r^M r^M 

obviously vanishes. Hence, the horizon area is left unaffected by the perturbation. 

B. The Time Dependent Problem 

When the sources are moved slowly from infinity into the immediate surroundings of the black hole horizon, we 
must retain the time derivatives in the scalar field equation 

^^A^^-m'%A-L'^ = ^. (4.5) 



(r - M)2 dry ' dr 

In analogy with the derivations in the previous sections, one should look for solution of the form 



$ = 5R/ ^ Q™(cc>)/,(a;,r)F,„(0,^)e-^'^*. (4.6) 



£=0 m=-i 



As usual, it is possible to transform the last equation with the aid of the transformations r* = r + 2M h\{r / M — 1) 
M/(r/M - 1) and Hi^{r*) = rfe{u;,r) into 

+Veir*))He^^u;^Hi^, 



dr 

mV /2M 2AP ^ £{£+!) 
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As in the non-extremal case, when the sources are moved through the potential barrier into the region where uP' is 
much greater than the potential, one can look for a solution of the form 

Hg^(r*) = exp(±icjr*) x [1 + 0(1 - Mjrf]. (4.8) 

Again, the Matzner boundary conditions select the sign in the exponent as negative. Hence, a typical term in <I> is as 
in Eq. ( 3.10D (with M equal to Q). As with the non-extreme case it can be easily verified that the invariant Va(f>V"<f> 



is bounded on the horizon, and is generally of order 0{uP). Finally the term TapN°-N^ reads 

Taf3N°'N^ = {d<^/dt + (1 - M/rYd<^/drf (4.9) 



and, after substitution of $ as in Eq. ( 3.10 )(with M equal to Q), it is plain to see that our former conclusion concerning 
the non-extreme case is recovered. 

V. DISTURBING A STATIC BLACK HOLE IN A DE SITTER UNIVERSE 

As a last example of processes in which the horizon area of a static black hole is left unchanged by an external 
scalar disturbance, we now consider the case in which scalar charges are distributed in a de Sitter universe with a 
Schwarzschild black hole embedded inside it. This spacetime can be described by the metric 

ds^ = -(1 - 2M/r - -r^)dt^ + [l - 2M/r - -r^Y^dr^ + r^{de^ + sinO^dip'^) (5.1) 
3 3 

where M is the mass of the black hole, and A is the cosmological constant. This spacetime is characterized by two 
horizons which correspond to the two positive roots of the cubic equation, gtt{r) = 0. If QAI^A < 1 there are three 
real roots. One can be attributed to the black hole horizon, rsH = fBH{M, A), while the other one is the cosmological 
horizon, rc = rc{M^K). The third root tq is always negative and thus unphysical. If 9M^A — 1, the two horizons 
coalesce and are located at r = (SAf/A)^/"^. In the case where QAf^A > 1, we have the unphysical situation in which 
the mass of the hole exceeds the energy associated with the cosmological constant. 

The shape of the effective potential that a scalar converging wave experience when it is propagating in this spacetime 

is 



y.(r*) = I 1 - ^ - ^r^ ) I ^ - ^A + 1^:-^ ) . (5.2) 



/2A/ 2 ^ ^(£+1) 
3 ' A r3 3 

In order to test our hypothesis regarding the adiabatic invariance of the black hole horizon area, we must consider the 
propagation of waves in the region between the black hole horizon and the location of the maximum of the potential 
barrier. This is due to an assumption that the squared frequency of the waves is small compared to the height of the 
potential barrier. Thus, waves which propagate outside this region will have no effect on the black hole's horizon. 
This problem was dealt with in px[ | and a generalization to the charged case was presented in the previous sections. 
Therefore, we shall concern ourselves only with the effect of scalar waves on the cosmological horizon rc, namely the 
effect of the waves in the region between the potential barrier and the cosmological horizon. 

For brevity we consider only the time dependent problem. As usual, we can solve the scalar field equation by a 
separation of variables. Thus, after expanding the time dependence in Fourier modes and the angular dependence in 
spherical harmonics, the radial part of the wave equation can be written in term of the difference of the roots of the 
cubic equation (?tt(r) = ^ = tq — tq, 6 = rc — rg^, and the radial location of the cosmological horizon compared 
to the radial distance to the cosmological horizon, namely e = — r, 

^rc-e)A— +^(^+1) + ^^-^ i? = (5.3) 



de \ de 

where A = A(7 — e)((5 — e)e/3. Note that 7,(5 and e are all positive semi definite by definition. 

In light of the previous derivation it is obvious that the problem has two competing scales; the proximity to the 
cosmological horizon, denoted here by e, and the distance between the two horizons, denoted here by (5. In the 
beginning we will assume that the black hole horizon is far from the cosmological horizon and that the sources are 
placed in a close distance from the cosmological horizon, namely ruH <^ r < rc- Later we will consider the possibility 
that the two horizons are on the verge of coalescence, (5 — > 0. 



Therefore, Eq. (5_^) can be approximated: 
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^(PR , dR ^ fZrcuj^'^ 



This is an Eiiler equation with the two independent solutions R{t) = N exp(±z77 In e), where TV is an arbitrary constant 
and rj — 3rcw/A7(5. Now, since the cosmological horizon is situated at e = 0"*", the solution with the plus sign should 
be discarded while the solution with the minus sign should be retained in accordance with the fact that only ingoing 
waves exist at the one-way membrane r = rc- Therefore, <i> consists of terms such as 

Fe(cos6') cos V', ijj = Lo{rih\€ + t) — raip. (5-5) 

It can be easily verified that the above solution does not cause any divergences in the invariants of the curvature. 
Hence we can turn now to the problem of verifying that the increase in the area of the cosmological horizon is at 
the same rate as the one found in the previous examples, namely proportional to o;^, while keeping in mind that the 



perturbation to the met ric is proportional only to lo (see Sec. III). To do so, we just have to calculate expressions 



such as the ones in Eq. ( 3.12 ) with N°' appropriate to the new spacetime. 



Tc^hN'^N'^ = (^^/dt+{l-2M/r-Y'^)d^/dT^ . (5.6) 



If one now substitutes a <I> made up of a single mode like in Eq. ( |5.5| ), one concludes that 

TapPl'^ cx w2p2 gjj^2 ^ ^_ . . . ^ (5 7) 

This result agrees with the ones obtained in the previous situations. 

Consider now a series of cases in which the separation between the black hole horizon and the cosmological horizon 5 
is gradually diminished. As long as (5 > 0, the potential barrier maintains its shape, and a weak scalar wave impinging 
on the black hole horizon can have no effect on the cosmological horizon and vice versa. If this is so, then our 
conclusion still holds: the black hole and cosmological horizons are left unchanged by the perturbation. In the limit 
where S the potential barrier, which separates the two horizons, disappears and the wave can either propagate 
between the innermost side of the horizon and the singularity, or in the region beyond the cosmological horizon. Both 
cases will be left for a future discussion. It is important to understand that S is not a dynamical parameter in the 
sense that for a given spacetime S is determined uniquely. However, it has a significant influence on the magnitude of 
the sources that may be used in the adiabatic process. This concludes the example. 



VI. KERR BLACK HOLE DISTURBED BY AN AXIALLY SYMMETRIC ARRAY OF SCALAR 

CHARGES 

In the previous examples we considered a static black hole disturbed by quasistatic perturbation. Now we shall 
study in brief the possibility that the hole itself is rotating. To simplify the calculation, we assume that the scalar 
charges are arranged in a configuration with cylindrical symmetry. We work in Boyer-Lindquist coordinates with the 
metric 

. , 2Mr, , Mar sin 9"^ 
ds^ = - 1 ^)dt^ - 4 dcjidt 

p2 J ^ A 



■ sind^ ( + + — sin 6*2 ) dtp^ + ^dr^ + p^dO^ (6.1) 



where A — r^ — 2M r + and p'^ — r'^ + cos 9^ . M is the mass of the hole and a is its specific angular momentum. 

As with the static spacetime case, we begin by writing down the scalar field equation in Kerr spacetime. The 
details of the solution of the scalar wave equation in a Kerr background without any assumptions on the scalar field 
can be found in [T^ . Assuming that the scalar field is time independent and has axial symmetry, one obtains an 
equation which is remarkably similar to the equation derived for the scalar field in the static spherical spacetime. The 
arguments we used in the static spherical spacetime are still applicable here. Therefore, we just write the physical 
solution 



$ - 3?^ QF, [:^^f=^ ] Pi{cos9). (6.2) 
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As before, it is possible to use the complex coefRcient Ci in the matching of the solution to any charge configuration 
with said symmetry. 

Now, in order to draw a conclusion regarding the adiabatic invariance of the horizon area, we look at the term 
T^N°'N/j where N" is proportional to the tangent of the null generator of the horizon which in Boycr-Lindquist 
coordinates reads [Q 

N* = 1, 
7V = ^ 
= 0, 

N'^ - (6.3) 

Since is nuU, namely N^'Nr = -{N^Nt + N'l'N^), we can write 

T^Nl^Nc, - (T/ - T^)N*Nt + {T^ - T^)mN^. (6.4) 

However, because of the symmetries, the terms T/ - T^: and - are both equal to -V^^V*. Thus, T^N^Na 

can be reduced to (N^d^/dr)^ — (d^/dr A/{r^ + a^))^. Since $ is analytic on the horizon, defined by A{r-H) — 0, 
the horizon area is left unchanged by the perturbation. 

The story becomes more interesting when one considers the time dependence of the sources. As usual, we make a 
separation ansatz which is given by the eigenfunctions appropriate for an axially symmetric and stationary background 
geometry, namely, ^e^i — Reuiir)Qe{6)e^'"^^ . By the usual argument we find the separated homogeneous equations, 
with separation constant p8[, 



dr \ dr 

-^ii (sine^)+aWcose'ee = ~XiQi. (6.5) 
sm u dO \ dO J 

The second equation is the same as the flat-space angular spheroidal equation with m = 0. Its eigenvalues are 
denoted by Xi and its eigenfunctions Se{~iauj,cos6). These eigenfunctions form a discrete set and go over into 
Legendre polynomials in the limit auj — s- 0. The integer £ has its standards range, but the eigenvalues Xi cannot be 
analytically expressed in terms of £ |l^] . 

The radial equation can be written as a one-dimensional equation with an effective potential by defining a new 
radial function and a new radial coordinate. 



He^ir*) = \/r^ + ami{u;,r), dr* = ——dr. (6.6) 



We then find 



dr 



*2 



T. / *N ^ 2 2. Sr^-AMr + a"^ 3r^A , ,^ . 



This potential depends nontrivially on the energy of the field. Now, from Eq. (6.7) it is obvious that in the 
proximity of the horizon, defined b y A (r-H) = 0, the potential is negligible with respect to the "energy" of the field 
Lo^. Accordingly, we can solve Eq. (|6.7D to first order in jFj/w^, and write the solution 



Hg^ = exp(±ti>r*) + (higher order terms) . (6-8) 

We choose appropriate boundary conditions for the solution of Eq. (|6.7]), analogously to Matzner's choice jl^ for the 
Schwarzschild case (a = 0). Thus, the minus sign in the w-mode expansion of $ and in Eq. (6.8) is singled out. This 
is a result of the demand that only ingoing waves exist at the one-way membrane r* —oo. 

Before undertaking the task of verifying that the perturbation does leave the horizon area unchanged, we have to 
verify that the solution is indeed physical, namely invariants of the curvature are bounded on the horizon. To do so 
we look, as before, at terms such cLS Vq$V'^$. a direct ca-lculation of the last term delivers the following result: 
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1 



c?ijp' sin il)^ sin d^Si^ + wr sin 2ipS£o + cos -0^ 



^0 



(6.9) 



where ip = uj{r* +t). Surely, the last result indicates that the solution we have found does not cause any singularities. 
Thus, the solution can be taken to be physical. 



Finally, following the chain of reasoning that led us to the derivation of Eq. (6^), we find that 



(6.10) 



This finding is in accordance with the results obtained in the previous sections. Thus we find yet another example of 
the adiabatic invariance of the black hole horizon area. 



VII. BLACK HOLE DISTURBED BY ELECTRIC CHARGES 



Here we shall examine the influence of an electromagnetic radiation on the black hole horizon area. In order to 
sidestep the problems of electromagnetic interactions between the black hole and the charges which produce the 
electromagnetic waves, we assume that the black hole is neutral. In the first subsection we formulate the problem. 
We make use of the possibility to describe the Schwarzschild space-time as a conformally flat spacetime, occupied by 
a "medium" which has a non-trivial permeabilities. In the following two last subsections we discuss the effect of the 
electromagnetic radiation on the horizon area in the static case and in the dynamic case. 



A. Formulation of the problem 

Following the footsteps of ||l^, Maxwell's equations can be cast into a simple form if we write the Schwarzschild 
metric in isotropic coordinates. Consider the coordinate transformation p defined hy r = M + p + M'^/4p and let 
f{p) = 1 — 2M/r{p) and u{p) = r{p)/p, then the metric reads 

ds^ = -f{p)dt^ + u{pf{dp^ + p^de^ + p^sinO^difi'^). (7.1) 

As far as the electromagnetic phenomena are concerned we can think of the space-time as Minkowskian, occupied by 
a "medium" which is characterized by dielectric and magnetic permeability tensors etk — Pik = n{p)Sik where 



n{p) EE u{p)/^J{p)- (7.2) 

n(p) can be viewed as an index of refraction induced by gravity. Note that n{p) is everywhere positive definite, and 
diverges as r ^ 2M. Maxwell's equations can then be written as 



± 



dt ' 

^ • {nP^) = 0, (7.3) 

where and H are the electric and magnetic fields in the medium respectively. The constitutive 

equations can be written as 

^ = n{p)'P, 

~§ = n{p)li. (7.4) 

Note that due to the static nature of the spacetime these constitutive equations do not include "gyroscopic" terms 
such as X H, where U represents the vorticity of the timelikc killing vector. V is taken in the last 

equation as the usual gradient operator in 3-dimensions, pd/dp + 0(1/ p) d/d9 + ip{l/ psin9) d/dip. p, 9 and ip are 
unit vectors in the radial, polar and azimuthal directions respectively. 
Consider now the representation 

^^{lf,t)^J2jd^^~'"' E F^^^\p,u.)y'':\p), (7.5) 

k a-—0,e,fi 
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k is an alias for (£,m). £ and m have the usual range. 1^^ and 'P*[^'' are the transverse and ^ the longitudinal 
vector spherical harmonics | pO| . They may be expressed by 

=py.^pYk{p). (7.6) 
When one substitutes this in Maxwell's equations, Eq. ([z^), one gets 

d 



dp 



d 
dp 



(pF^^^^) =±u:pn(p)F^^^\ 
_ ./KITT)F^^'^ = ±u;pn{p)F^^''\ 



An indication of the relevance of the solution to the physical problem at hand can be obtained by examining the 
invariants that may be constructed from the energy-momentum tensor of the perturbation 

= 4^ {F,<^nf'"6^p - \s;Pcpn-') , (7.8) 

where we define the usual field strength, Fap and the pseudo-tensor Ti.°'^^ = ^/—gF"^ {—D,7i). Note 

that \/—g = /n. Every invariant of this kind must be bounded for a physically acceptable solution. The analysis 
can be simplified much using the following arguments. may be seen as a 4 x 4 real matrix [2l| ]. Thus, any 
invariant constructed from is actually the trace of some power of T. For example, = TrT = 0, T^T/^ — TrT^, 
T^Tl^T^ = Tr etc. Now, in this matrix notation can be written as 

^ M - ^TrM 1 ) (7.9) 



47r^/^ V 4 

where = F^TC^ or in matrix notation M = F • ?i and 1 is a 4 x 4 identity matrix. The trace of any power of T is 
simply the trace of some polynomial function of M, where the polynomials coefficients may be powers of Tr M. Now, 
since the trace operator is linear, it is sufficient to look at expressions like TrM, TrM^, TrM'^ etc. These in turn are 
simply the sum of the eigenvalues of M raised to the appropriate power. Hence, a prerequisite for the boundness of 
the invariants of the curvature is that the eigenvalues of M vanishe at least as fast as 1/n as we approach the horizon. 
It turns out that M has only two distinct eigenvalues which have each a two-fold degeneracy. These are 

Ai = J + 0'\\)\ \2 = \ (ll i^^ll (7.10) 



where . Now if one uses the definition of the vector spherical harmonic Eq. ([7^) it is plain 



to see that i^l^' • ^1-2^ ~ "^fei^ ' "^ia^ for every k, and since the vector spherical harmonics form an orthogonal base, 
II II can be reduced to 



1+1^2)* 



lli^^ll' = ^ H / dLUidLU2e-''-^'- 

fci ,k2 

f n,n.F,t(°V,t(°) + Vn, . Vn^ ^ Ft^^^^F^'] ■ (7.11) 

\ tT=e,p, ) 

In the following sections we will use this expression to verify that our solution is physically meaningful. 
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After the solution is shown to be physical, we will test whether the horizon area is left unchanged under the 
disturbance. This will be confirmed if the term AirTajjl^l^ vanishes on the horizon. In Sec. Ill P, the null generator 
of the horizon was shown to be proportional to the vector N"-. Now, under the coordinates transformation to the 
isotropic coordinates system, the vector iV" transforms to 

iV" EE (1, 1 - 2M/r, 0, 0) ^ (1, 0, 0). (7.12) 

Thus, a straightforward calculation leads to the following result: 

^TTTa^prf oc + iF^+) {Fg- - iF^-) (7.13) 

and after substitution of from Eq. (|7.5|), one obtains 

fe ^ 

where is the differential operator d/d9± {i/ sin 6)3/ dip. All the above is general. In the following subsection we 
shall discuss the static problem. Finally, in the last subsection we shall retain the time dependence of the fields and 
solve the dynamical problem. 



B. The Static Problem 



When the electromagnetic fields produced by the charges are static. Maxwell's equations [Eq. ([7.7D] can be reduced 



to 



yMiTT)F^''^^ 



,±(0) 



dF, 



±(0) 



im 



(ln(p2„))Ff^°^+^(^+l)F, 



±(e) 
k 



0, 
0, 

0, 
0. 



dp dp 

The first two equations simply state that F^'"^^ = for every k. The last two equations can be combined to give 



(7.15) 



F, 



±(0) 



(7.16) 



i^^^"^ is given by the second equation of Eq. ( 7.16| ) after the first equation is solved for F^^'^K 

After the substitution of 'n(p) and the introduction of a new function A'^'^'^^ = pF^^'^\ we can write the equation 
for F^^'^'^ in the following form 



£ u 2 



^ + Wi + i)f- ^ ^ 

(p+fKp-f) dp +^'^' + '^^ p 



(7.17) 



Subsequently we will be interested in the behavior of th e electromagnetic fields in the proximity of the horizon. 
Therefore, we look for an approximate solution to Eq. ( 7.17j ) near p = M/2. We introduce a new variable x = p — M/2 

and a new function X^^'^\x) = xz^^'^^x). Then Eq. (7.17) assumes the following form 



d'^z^^"^ 
x^ — A ha; 



dz 



dx'^ 



dx 



1) 



(7.18) 
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with a = 2{£{^-\- 1))^^"^ /M. This equation is the famihar Bessel equation of or der 1 with the two independent 
solutions Ji{^/ax) and Yl[^/ax). Convergence of the curvature invariant, Eq (7.11), forbids us from considering the 
Yi{^/ax) solution as a physical one since it diverges at x = 0, namely at the horizon. As in the former cases, this 
divergence is unacceptable since it wiU induce divergences in the invariants of the curvature, thus disqualifying the 
perturbation method. 

Using Eq. (7.16), it is now possible to show that 



F, 



±(m) _ 

±(e 



0, 

X 



k - ^C^MVax), 



[Jo{y/ax) — J2{\/ax)) 



(7.19) 



Hence, according to Eq. ( 7.13 ) and Eq. ( 7.14 ) (with uj 



behavior of Fj^^"' near the horizon. Now, for a; ^ 0, Ji{y/ax) is of 0{x). Therefore, i^^'"^"* is of 0{x 



0) the behavior of inTapl^l^ is governed solely by the 

2 



and F^'"^^ is 



of 0{x). This implies that AnTa/sl^l^ is of order 0{x'^). Thus the horizon area does not change under the action of 
the static electromagnetic fields. 

Since the adiabatic invariance of an object is primarily meaningful in the context of dynamic processes, we turn 
now to deal with the dynamic problem. Basically, we will consider the process in which the charges are brought 
adiabatically from infinity through the potential barrier of the black hole into the immediate surroundings of the 
horizon. This will be dealt with in the next subsection. 



C. Moving the Charges Around 

First we introduce a new coordinate r* defined by dr* /dp = n{p). Then r*{p) can be identified with Wheeler's 
"tortoise" coordinate |Q; r* — r + 2Af In (r/2A/ — 1). Moreover, as usual, we define new radial functions 

X^'^''\r* ,uj) — pF^^'^K In terms of these new radial functions the radial parts of Maxwell's equations reads 



dr* 



±(e) _ , .dX,^ 



±(m) 



dr* 



A±(") = T(i/^)^I(7TT)^. (7.20) 

We are interested in a solution in the region where lj^ is much larger than the potential barrier, namely in the 
region where 

.^»ri-^)^. (7.21) 



r 



An approximate solution is 



A±('')^C±(a;)e 



A±(°) ^ ^^^(TTT)^^^^^^^^. (7.22) 

U! pn 

A useful by-product of the fact that 1 /n vanishes on the horizon, is that the longitudinal component of the electro- 



magnetic fields vanishes as y^l — 2M/r when r — > 2M. 



After substitution of A^'^^'^^ from Eq. ( 7.22 ) in Eq. ( 7.11 ), it is evident that the summation over a vanishes identically, 
regardless of the magnitude of the frequency dependent coefficients, C^(cj). And as mentioned, the longitudinal term 

F^'-^'f^'"'' is proportional to 1/n^. Hence, the multiplication of this term by n does not cause any trouble. In fact as 
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we approach the horizon the contribution from this term to the perturbation of the metric becomes neghgible. This 
shows that every invariant of the type we considered vanishes on the horizon in spite of the divergence of n there. Of 
course, globally these invariants may still have an effect on the geometry. In particular the geometry of the horizon 
can get disturbed. 

Now, since the perturbation is supposed to be weak, we are limited to small frequencies, namely long wavelengths 
compared to the size of the hole, A >> 2M. Additionally the coefficients C^{uj) must be bounded everywhere and are 
assumed to be smoot h fun ctions of co. A s before, consider a solution composed of a single mode. Then, substitution 
of this solution, Eq. ( |7.22| ), in Eq. ( |7.1l| ) gives 



£{£- 



cos 



= o 



n 2N 



(7.23) 



where ip = uj(r* +t) — mtp. It is evident from above that at a small radial distance from the horizon the electromagnetic 
fields induce an oscillatory disturbance of 0([Cfc(ci;)/a;]^). The disturbance vanishes on the horizon. The limit w — s- 
was dealt with in the previous subsection and was found to exist. Now, how does the horizon area change when the 
weak perturbati on is turned on ? Co nsider again a solution composed of a single mode. Then, the substitution of this 
solution in Eq. ( 7.14 ) and Eq. ( 7.13| ) gives 



oc 



2 cosip 



fdPi 

\ de 




(7.24) 



Obviously, this result states that the perturbation to the horizon area is a quantity of higher order in w than the 
perturbation to the geometry as quantified by the invariants calculated previously. Thus horizon area is again seen 
to be an adiabatic invariant. 
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